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Abstract 

The heat current across a quantum harmonic system connected to reservoirs at different temper- 
atures is given by the Landauer formula, in terms of an integral over phonon frequencies a;, of the 
energy transmittance T{oj). There are several different ways to derive this formula, for example 
using the Keldysh approach or the Langevin equation approach. The energy transmittance T{oj) 
is usually expressed in terms of nonequilibrium phonon Green's function and it is expected that it 
is related to the transmission coefficient t{uj) of plane waves across the system. In this paper, for 
a one-dimensional set-up of a finite harmonic chain connected to reservoirs which are also semi- 
infinite harmonic chains, we present a simple and direct demonstration of the relation between 
T(w) and t{uj). Our approach is easily extendable to the case where both system and reservoirs 
are in higher dimensions and have arbitrary geometries, in which case the meaning of r and its 
relation to T are more non-trivial. 

PACS numbers: 65.40.Gr,05.40.-a,05.70.Ln,44.10.-|-i 



I. INTRODUCTION 

The Landauer formula gives an exact expression for the current (energy and/or particle) 
in "non-interacting" quantum systems coupled to reservoirs kept at different temperatures ( 
and/or different chemical potentials ). By "non- interacting" one refers to systems described 
by quadratic Hamiltonians. It thus includes harmonic crystals where one considers energy 
transport by phonons, and tight-binding Hamiltonians where there is transport of both 
charge and energy by electrons. The formula for phonon heat current across a harmonic 
crystal connected to heat baths at temperatures Tl,Tr is given by 



where the quantity T{u)), which we shall refer to as the energy transmittance, can be ex- 
pressed in terms of appropriate "nonequilibrium" phonon's Green's functions and f{u,T) = 
-^jl^^hw/ksT _ ^j^g thermal phonon distribution function. This Landauer formula for 

phononic heat current has been derived rigorously using the quantum Langevin equation 
approach [U |2] as well as the nonequilibrium Green's function (NEGF) approach [3l H]. 

Landauer's original idea was to think of conductance in terms of transmission or scat- 
tering of plane waves and for the case where the reservoirs or "leads" are one dimensional, 
it is expected that T{(jj) is related to the transmission coefficient t{uj) of plane waves [5]. 
For the case of electron transmission, T(w) can again be expressed in terms of nonequi- 
librium Green's functions [6] and the relation to the transmission coefficient was directly 
demonstrated through the work of Todorov et al using scattering theory [7J . For the case of 
phonons we are not aware of an explicit proof of this relation and that is the main objective 
of this paper. Here we consider a general one-dimensional finite harmonic chain coupled 
to reservoirs which are themselves semi-infinite ordered harmonic chains and give a fully 
quantum-mechanical derivation of the relation between T and r and show how the NEGF 
current formula can be obtained from the transmission coefficient. 

In our derivation we first note that in the NEGF approach the energy transmission T is 
expressed in terms of a Green's function. This Green's function can be expressed explicitly 
in terms of a product of 2 x 2 matrices. On the other hand the transmission coefficient 
can be computed by constructing appropriate scattering states and this can be done in two 
ways — (i) a direct solution of the discrete wave equation which again gives r expressed 
in the form of a product of matrices or (ii) by using the Lippmann-Schwinger scattering 
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theory to evolve reservoir normal modes and this gives r directly in terms of the Green's 
function. From the forms of these expressions we directly obtain the required relations. 
We note that for the case where the reservoirs are not one-dimensional chains, but have 
arbitrary geometries pHTO]. the NEGF expression for T{uj) still has the same form but it is 
not clear as to how one should compute r and how exactly it is related to T. In this case 
the approach using Lippmann-Schwinger scattering theory can still be used to arrive at the 
required relation. A model similar to ours was studied recently by Zhang et al [11] in the 
context of interfacial thermal transport in atomic junctions and the relation between the 
NEGF formula for energy transmittance and the transmission coefficient was established 
numerically and also exactly for the special case of a single interface. 

The plan of the paper is as follows. In sec. |TT] we first define the model and state some 
general results for the heat current given by the formalism of nonequilibrium Green's func- 
tions. We then give an explicit expression for the form of the Green's function appearing 
in the energy transmission formula. In sec. |III| we consider the transmission of plane waves 
across the system and, using two different approaches, obtain the form of the transmission 
coefficient. The transmission coefficient can also be expressed in terms of the same nonequi- 
librium Green's function and using this we write the relation between it and T{(jj). This 
relation is then used in sec. |IV| to derive the Landauer formula for heat current. Finally we 
discuss our results in sec. |Vl 

II. ONE-DIMENSIONAL CHAIN CONNECTED TO ONE-DIMENSIONAL 
BATHS 

Consider the set-up in Fig. ([T| where a one-dimenional (ID) harmonic chain with arbitrary 
spring constants and masses is connected to leads which are themselves ordered harmonic 
chains. Special cases of this setup have been discussed earlier by various authors [TTHT6] in 
the context of heat conduction. Let us assume that the system has Cartesian positional 
degrees of freedom {x;}, / = 1, 2 . . . , with corresponding momenta {pi}- These satisfy the 
usual commutation relations = i^Si,m and = [pi,Pm] = 0. Similarly the left 

reservoir degrees of freedom are denoted by {xa,Pa}, a = 1, . . . , and the right reservoirs 
by {xa',Pa'}, a' = I, . . . , Nfi. We consider our system plus reservoir to be described by the 
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FIG. 1: Schematic of the set-up considered in the paper. The system consists of a harmonic 
chain of N particles for which both the particle masses and the inter-particle spring constants take 
arbitrary values. The system is sandwiched between two reservoirs which are ordered ID harmonic 
chains with different mass densities and elasticities. The coupling constant between left reservoir 
and system is ko and between right reservoir and system is /ctv- 

full Hamiltonian 

1=1 ' 1=1 

^ 2 2 

a=l ^ a=l 

pI' I knjxc,' - Xg'+i)'^ fcjv(x„/=i - xnY , . 

+ 2^2mn^ 2 + 2 ' ^"^^ 

a'=l " a'=l 

where we assume Xa=NL+i = Xa'=NR+i = 0. The system masses {nii} and spring constants 
{ki} are assumed to be arbitrary. The left (right) reservoir particle masses are all taken to 
be {mfj) and the inter-particle spring constants are taken to be fc^ (/c_r). To ensure a 
uniqe steady state We will always assume that the reservoirs are chosen to have sufficiently 
broad bandwidths compared to the spectrum of the system [H [18] . The above Hamiltonian 
can be written in the canonical form: 

"H = Tis + T^L + 'Hr + T-Lls + T^RS , (3) 



4 



where 

17 , sr^ ki{xi - xi+if kpxi^ k^XN"^ 

- l^^ + l^ 2 + — + ' 

1=1 1=1 

U _ \^ pI , kiiXa-Xa+iY koxl^^ 
a=l 

u _ pI' , knix^, - Xa'+if , k^xl,^^ 

a'=l 

T-LlS = —koXa=lXi, T-Lrs = —kNXa'=lXN ■ (4) 

Using the vector notation Xj = {xi,X2, ■ ■ ■ ,Xn), PJ = {PiyP2, ■ ■ ■ ,Pn) and similarly 
Xl, Xpi^, Pl, Pji, the different parts in the above Hamiltonian can be written as 

ns = ^Pj M,' Ps + ^Xj Ks Xs , 
Hl = \PJ M^' Pl + ^Xj KlXl , 

'Hls = KsL Xl , Hrs = Xg KsR Xr , 

where Ms, Ml, Mr and Ks, Kl, Kr denote respectively the mass matrix and the force- 
constant matrix of the system, left reservoir and right reservoir, while Ksl and Ksr denote 
the linear coupling coefficients between the two reservoirs and the system. In our case Ksl 
is a X X Nl matrix whose only non-zero element is [1^5/^]i i = kp, while Ksr is a X x 
matrix whose only non-zero element is [Ksr\n,i = k^- 

Expression for steady state heat current: We now consider the situation where 
at some distant past time {t < to) the two reservoirs are uncoupled from the system and 
are separately in equilibrium (and described by canonical distributions) at temperatures Tl 
and Tr respectively. At time to we start evolving the system plus reservoirs with the full 
Hamiltonian in Eq. ([s]). Eventually we set the reservoir sizes NL,Nji — )■ oo and to ~^ — oo. 
The system reaches a nonequilibrium steady state at finite t. Note that we have included 
terms involving the coupling coefficients fco, /cAf in the isolated reservoir Hamiltonians. As has 
been discussed using various approaches [H El E] , the steady state current can be expressed 
using the following phonon Green's function: 

^ -Msu- + Ks-^t-J:%' 
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where the self-energies S^, can be expressed in terms of the isolated reservoir Green func- 



tions g^(a;) = [ -Mi(a;±ze)2 + Ki]-i 



' Br 



[—Mji{uj±iey + Kji\ ^ and the coupling 



matrices Ksl, Ksr- The self energies are given by H^iu) = Ksl sti^) -^sl^ — 
KsR g^ico) . Defining TLico) = Im[ S+ ] , r^(a;) = Im[ S+ ], we find P El S] that 

the steady state current is given by the formula in Eq. ([T]) with 



Tito) = ATr[G+iLo)TL{co)G.{Lo)T Rico)] . 



(6) 



For our one-dimensional system, we note that G^,S^,S^ are all x matrices. The 
only non-zero elements of and are respectively [S^]i^i = A;o^[s'^]i,i =: and 



/cAr^[gi^]i,i =: S^. Let us define Tl = Tr = /m[S^] . Hence the 



expression of T reduces to: 



(7) 



with the matrix 



r-l 



where Z = —Ms^"^ + -^s ~ ~ ^ tri-diagonal matrix. 



The bandwidth of the two baths are different {2^JkL/mL and 2^JkR/mR for the left and 
right baths respectively), so the conduction of heat across the system will have contribution 
only from the overlapping part of the bandwidths. 

Explicit forms for G^ i^,T r. Using methods described in [T7j we now show that 
the Green's function element occuring in Eq. ([7|) can be expressed in terms of a product of 
2x2 matrices. We also obtain the explicit forms of F^,, Tr for our particular model. The 
matrix Z has the form 



Ziuj] 



ai - T^Ku) -ki ■ ■ ■ 
—ki a2 —k2 












where ai 



■ ■ ■ —kN-2 ClN-l —kN-1 

• • • —kN-i cln — Sr(i^) 

ki + ki^i - miu^ , l = l,...,N. 



(8) 



Taking the inverse of this matrix, we get 



N-1 

n k 

1=1 



Ai. 



N 



(9) 
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where Ai^n is defined as the determinant of the matrix Z{uj). Let us also define Di j as 
the determinant of the sub-matrix starting with the l-th row and column and ending with 
the j-th row and column of the matrix —MsOj"^ + Ks- From the tri-diagonal form of the 
matrices, it is easily shown that 



Ai 



N 



ai - S J ) [ (ajv - S^) D2,N-i - kN-iD2,N-2 ] - ki [ (aat - 'L']^)D?,^n-i - kN-iD^^N-^ 



The elements Di j satisfy the recursion relation Di^j^ 
h,N-i 





1 







aiDi+i^N-kfDi+2,N for / = 1, ...A^-2, 



and Di,N_^ = aiDi+i^^^i — kfDi^2,N-i ioi I = 1, N — 3. In matrix form these give 



which holds for / 



h 



ai/ki 

i/h 



-ki 




l+l,N 



l+2,N 



-Di 



+2,N-l 



(11) 



l,..,iV 



D 



N,N-l 



1, D 



N+2,N 



3. Using these relations and further defining -Dat+i^at 
Dn+i,n-i = 0, we arrive at the result 



1=1 



where T 

Hence using Eqs. ( [I0p2p| we get 

N 

Ai,^ = (n^'; 



D2,N —D2,N-l 
N 



N 



1=1 

ai/ki -ki 
l/ki 



1 

l/k 



N 



(12) 
(13) 



1 -^t 



T 



1=1 



1 




1 


i/kl 







(14) 



We next find the explicit forms of S^, Sjj, for which we need to evaluate the reservoir Green's 
function elements [g/^]i,i and [g/j]i,i- Consider the left reservoir. For the case ko = ki, it is 
simple to find all normal modes and hence compute the Green's function corresponding to 
the force matrix Kl = (say). One gets 

1 



-Ml (co + ie) 



Hence [g^'^ 



Lm 



where fl„ 



mLir Jo 

2kL 



dq 



sin(g/) sin(gm) 

-{uj + itY + ni 



rriL 



[1 - cos(g)] . 



(10) 



We need the (1, 1)**^ element and the above integral gives [g/^^]i,i = e^'^/fci, where q is to be 
obtained from uj"^ = {2kL/mL){l — cosg). For the general case ^ ki, the Green's function 
can be calculated as follows. We write Kl = + AKl where AKl is a perturbartion 
matrix whose the only non-zero element is AKn = k^ — k^. From the definition of the 
Green's function g;^ = [-Ml{u} + ie^ + Kl + AKl]''^ we get 



Taking the (1, 1)* element of the above equation gives 



1,1 



[g^] 
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3«g 



1 + [g°+]n [AKl]u kL + (/Co - A:L)e*'^ 



and similarly 



="7 



1,1 



kR + {kN - kR)eii' ■ 
Using the definitions given earlier we derive the following expressions: 



(15) 



(16) 



(17) 



kle^'i 



^2 iq' 



^ kL + iko-kL)e'i' ^ kR + {kM-kR)ei<i' 



kQkLsin{q) 



Ikn-kL + kL e-'^^ 



2 ' 



r 



R 



kjfkRsin{q') 



\kN -kR + kR e-^i' 



2 ' 



(18) 



where q,q' are respectively obtained from the relations u"^ = (2/ci/mL)(l — cosg) = 
(2A;^/m/j)(l — cosg') and rL,rR are non-zero only when both g, g' are real. Plugging in 



the expressions of S^^ and in Eq. (|14|), we obtain from Eq. (|9 



1 



1 




f 


kN 







1 


fc£ + (fco-feL)e*9_ 





1/kN 

















(19) 



III. SCATTERING STATES AND TRANSMISSION COEFFICIENT 



For our model the equations of motion correspond to the discrete wave equation for 
which we can construct scattering wave solutions. We will now construct solutions that 
correspond to plane waves incident on the system from the reservoirs. From these solutions 
we will obtain the transmission coefficient. In the following we will only consider the "right- 
moving states" which correspond to waves that are incident from the left reservoir. The 
"left-moving states" can be similarly obtained. 
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Let us consider a chain described by the Hamiltonian of Eq. ([2]) with an infinite number of 
particles in both the reservoirs. The particle displacements in the chain satisfy the equations 
of motion 



mixi = -{ki_i + ki) xi + ki_i Xi^i + k Xi+i , (20) 

where / = 1, . . . , refers to particles of the system, / < refers to particles of the left 
reservoir {i.e a > 1, as in the notation of Eq. [2]), and / > + 1 refers to particles in the 
right reservoir [i.e a' > 1). We note that these equations are valid both for the quantum 
representation, where the variables are Heisenberg operators, and also for the classical case. 
Corresponding to the above equations let us construct classical wave solutions ipi satisfying 
the equations 

mitpi = + ki) ipi + ki_i ^/„i + ki ipi+i . (21) 

In the left and right reservoirs these equations take the form of the discrete wave equations 

= ikL/mL){^a+i - 2V'a + ^a-i) for a > 1 , 

V'a' = {kR/mji){iJa'+l - '2^a' + ^a'-l) for > 1 . 

These have the following plane-wave solutions, 

"^"^^^ " (27rm )V2 """"^ ^ + ) for a > 1 , (22) 

where the wave-vectors q,q' G (0,7r) satisfy the dispersion relations 

tu^ = {2kL/mL) (1 - cosg) = {2kR/mR){l - cosg') , (24) 

and the normalization is chosen such that for r = {i.e no transmission) the following 
condition is satisfied: 

/■TT 

dq mLi)a{q)i)u{q) = (25) 





for any two points a, z/ on the left bath. The solution in Eq. (22) corresponds to a plane 
wave of wave vector g, frequency u that is incident on the system from the left side, part 
of this is then refiected with amplitude r, and a part transmitted across the system with 



amplitude r. We shall refer to r as the transmission coefficient and will now proceed to the 
calculation of this. As is well-known in quantum mechanics and wave-theory, the required 
scattering states can be constructed either by direct solution of the equations of motion in 



Eq. (21) or by the Lippmann-Schwinger scattering theory approach. We now present both 



these methods. 



A. Transmission coefficient from direct solution of the the wave equation 



For points on the reservoirs the plane wave solution has the form in Eq. (22). For 
< / < A^ + 1, let us write '4'i{q) = ■s;e~*'^*, Z = 0, .., A^ + 1 where the amplitudes si satisfy 
the equations 

niiu'^si = {ki_i + ki) si - ki_i si_i - k si+i , (26) 
and it is to be understood that / = 0, / = — 1 refer to a = 1, a = 2 respectively while 



/ = A^+1, / = A^ + 2 refer to a' 


= 1, a' = 


2. Hence we 


get the following recursion relation: 




k-lSl-i 


= fi 


kiSi 


(27) 




Sl 









where Ti is defined in Eq. (13). Using this recursively gives 



k-is-i 
So 



TaTT^ 



N+l 



kN+l^N+l 
Sn+2 



(28) 



We note that 



k-is^i 


1 


kL (e~2*« + re2^9) 




kN+lSN+1 


1 


kRc'i' 


So 


(27rmi)i/2 


(e-*g _^ re'") 




Sn+2 


(27rmi)i/2 






To = 


kL + ko- ■rriLUi'^ 
1 





, Tjy+l = 


kR + kN - mRU^ -k\ 
1 



T 



and hence Eq. ( 28 ) gives 



Tn TT-, 



N+l 



kRC^'i' 



r 



+ 



r . 



-iq 



,iq 



(29) 
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To solve for r we multiply the above equation by the row vector (1 — ^Le*'), and this gives 

—2ikL sin(g) 



T 



- kie^'i) To f f, 



kjie 



iq 



N+l 



,2iq' 



After some simplifications and using the form of G^j^ given in Eq. (19) we obtain: 

2ikLsm{q) kok^e-'^t+'i'^ 



-TV 

n' \ J-T 



(30) 



(/Co -kL + kLe-ii){kN -kR + kne-^i' , 
Now using the expressions for Vl^Vr in Eq. (18) and comparing with the formula in Eq. ([T]) 
we immediately see that the energy transmittance T{uj) and the transmission coefficient 
riu) are related as 



T{uj) = \t{uj] 



^ kn sm{q') 
ki sin(g) 



(31) 



B. Transmission coefficient from Lippmann-Schwinger scattering approach 

The Lippmann-Schwinger scattering theory approach in quantum mechanics starts by 
breaking up the Hamiltonian of a system into an unperturbed part and a perturbation. One 
then writes an exact scattering solution of the unperturbed part of the Hamiltonian and 
then uses this to obtain a solution of the full problem in terms of the perturbation and 
appropriate Green's functions. Here, using the notation of Eq. (|3]), we treat l-Ls + + I-Lr 
as the unperturbed Hamiltonian and 1-Lls + T^rs as the perturbation. 

Lippmann Schwinger theory: Let us use the notation \ip{q)) to denote the state for 
the wave-function il)i{q) satisfying the wave equation 



Ma;2 |^(^)) ^ K |^(g)) 



(32) 



where M and K are the mass matrix and force matrix respectively of the full chain (including 
system and reservoirs). Using the partition of the chain into the reservoir and system parts, 
these matrices have the following block structures: 



M 



Ms 
Mi 
Mr 



K 



Ks KsL KsR 
Kl 
Kn 



^SR 



(33) 
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Breaking K into unperturbed and perturbed parts we have: 



K = Ko + Ki 







Ks 













KsL 


KsR 


where 


Ko = 





Kl 

























Kr 













(34) 



Treating Ki as a perturbation we then obtain the following scattering solution of Eq. (32): 



1^) = \r) - Q K, m , 

where g+ = [-M{uj + tef + K]-^ 



(35) 



is the Green's function for the full chain and {ip^) is a scattering solution of the unperturbed 
system satisfying the equation 



(36) 



Construction of initial state: Let us first construct the right-moving scattering states. 
For this we consider the particular initial state where the left reservoir is in a normal 
mode with frequency u while the system and right reservoir degrees of freedom are at rest. 
Thus we choose ipfiq) = for / > and ipi{q) = ip^ ^ — 0; with a = 1 — I and 
satisfying the equation 

m^u;2^^(g) = J][K^]„,^^J(g) . (37) 

The form of can be read from Eq. (|4]) , and we then get 
niLOO^^M) = kd 2^^{q) - ^^^,{q) - ^tii<l) ] + '^a,i(^o - kL)^i{q) , a = 1,2,... ,(38) 

with the boundary condition V'o (o') = 0- For = the normal modes are given by ip^i^) ~ 
2ismqa/{27imLY^'^, where the normalization is chosen such that Jq dq fnLip^*a{q)ip^ (q) = 



Sa,u- For ko ^ k^ we can find the normal modes by treating the last term in Eq. (38) as a 
perturbation. We will require only ip^^^^q). The Lippmann-Schwinger approach is applied 
again, giving 



,r , , 2i sin(g) , n r +n i 2i sing , ^ ,^ 

^-l(^) = ,1/2 [1 - (fco - kL) [g+]ij = - ^^^^^y„ kLC [gt],,v 



(27rmi)i/2 



(39) 
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where the result in Eq. (16) has been used. Note that our choice of ipo imphes an incident 

wave e*'?"/(27rmi)i/2^ 

Scattering state: Since we want to finally find r, it is sufficient to compute the scattering 



(40) 



wave function only on the right reservoir. From Eq. (35) we get 

We now express the Green's function element —Q"^, in terms of the Green's function G"*" 
defined earlier in Eq. (|5|. We first write in a block-matrix form, with the different blocks 
representing the system and reservoirs. This matrix satisfies the following relation: 

^SR 



Ks KsL KsR 

-M^ {uj + iey + Kl 

-M^ (uj + ieY + K 



R 







GsL 


GsR 




I 








X 






G^R 







/ 







Grs 


G\l 


G^, 










/ 



(41) 



From this equation the following relations can be shown to hold [T]: 

1 



Gj(u;) 



-u'^MS + Ks 



Grs - SrKsrG 



This then gives us 

^l',l=l = [GRs]a',l=l = [gflla'.l^AfG^,! • 

Using Eqs. ( |39p4l ) in Eq. (|40]) we finally get: 

ipa'iq') = 2ifcofc7vA;Lsinge-^^[g^]i,i[g^]^/,iG+^/(27rmL)i/^ 



(42) 
(43) 

(44) 

(45) 



Now taking the {a', element of Eq. (15), with L replaced by R, and using Eq. (17) we 
get 



'■>r\»',i 



kR + (kN - kR)e^<i' 



(46) 



Using the explicit form of [g/^]i,i from Eq. (16) we finally arrive at the expected form of the 



transmitted wave function in the right reservoir 

^a'(9')=re^''°7(2^"^L)'/' 
13 



with r precisely given by the same expression Eq. (30) obtained in the previous sub-section 
by the direct solution of the wave equation. 



IV. EXPRESSION FOR THE ENERGY CURRENT IN EACH MODE AND A 
DERIVATION OF THE LANDAUER FORMULA 

We now use the definition of the heat current operator and show how it can be used to 
express the current contribution of each of the modes in terms of the transmission coefficient 
and hence the energy transmittance. This will lead us to a derivation of the Landauer 
formula. In the steady state the current is constant everywhere and we will evaluate it 
on the right reservoir. Between sites a' and a' + 1 the left-right current is given by the 
expectation value [2J ji = ( \kfi{va' + — Xa'+i) ) where we compute the average 

using the q^^ right-moving state obtained in the previous section. It is easiest to obtain this 
using second-quantized notation. The set of right moving and left moving states form a 
complete set. Denoting the left-movers by i'a'iq) with q G (— vr, 0) we note that they satisfy 
the completeness relation J^^dq(ma'Y^'^{mi,'Y/'^ip'^,(q)^,y/{q) = 5a' y- The displacement and 
velocity operators at the lattice sites of the right bath can be expressed in terms of the 
creation and annihilation operators 

h \^ 



x„/ = / dq[ — ] (agipa'iq) + alipa'*{q)) , 

aq-ipa'iq) - alilJa*{q)) ■ 



1 



Va' = - 

The operators a^j, aj^ satisfy the commutation relations [a^^, aj^] = S{qi — ^2) and, using the 
completeness relation, it can be verified that this ensures the usual commutation relations 
for the position and momentum operators. Using the above we get for the expectation value 
of the current for a right moving state: 

JLRiq) = ihkR{{alaq) + ^)['^a'+i*{q)^a'{q) - '^a'+i{q)'ipa*{q)] 

= \r\ [/K,T.) + -], (47) 

where in the last step we have used the form V'o'(Q') — tc'^"'''^' /{2TrmLY^'^, and the initial 
occupation probability of the state q is given by the left bath thermal distribution {o)qaq) = 
^^nwq/kgTL _ — f(^u^Ti). The total current transmitted from the left bath to the right 

14 



bath, is obtained by integrating over all q. After making a change of variables from q to 



uj = liki^jmi)^!'^ ^\ — cos(g) we get 

/■TT 

Jlr = / dqJiRiq) 
Jo 



-/ du hw- : r [f{uj,TL) + -] 

TT Jo kL sm g 2 

. p2{kL/mL)^/'^ ^ 

- dcohwnco)[f{co,TL) + -], (48) 



where in the last step we used Eq. (31) . From symmetry, the current flowing from the right 
bath to the left bath will be given by 



Jrl = - du huT{uj) [f{u,TR) + 



2' 

Hence finally we get for the net current: 

J=- dionwT{io)[f{uj,TL) - f{io,TR)] , (49) 

Jo 

where Um = niin[2(fci/m/^)-'^/^, 2(A;^/m^)^/^]. Observing that T is a symmetric function of 



Lu and vansishes outside the range u G (O,^^), we can see that Eq. (49) is equivalent to the 
Landauer formula Eq. Q . 



V. DISCUSSION 



In summary we have studied heat conduction across a ID quantum-mechanical harmonic 
chain, with arbitrary distribution of masses and inter-particle spring constants, that is con- 
nected to two other ordered ID harmonic crystals which have different mass densities and 
elastic constants. For this model we use two different approaches to demonstrate the relation 
T{ijj) = {kRsinq') |r((X')p/(/ci sing) between the energy transmittance T(w), which occurs 
in the Landauer formula for heat current, and the transmission coefficient t{u) related to 
passage of plane waves across the system. In the first approach we use the fact that the 
Green's function occuring in the expression for T has a simple representation in terms of 
product of 2 X 2 matrices. The plane wave solutions are then obtained by directly solving 
the equations of motion and a representation of r is obtained, again in terms of the product 
of matrices. The connection between T and r is then directly obtained. This approach can 
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be extended to the case of regular lattices using, for example, the techniques used in [19] for 
the representation of the Green's functions using matrix products. 

In the second approach it is not necessary to find the explicit form of the Green's function. 
One notes that the required plane wave scattering states can be obtained by using the 
Lippmann-Schwinger approach to evolve initial states which are eigenmodes of either one 
of the reservoirs and are initially localized within the reservoirs. The Lippmann-Schwinger 
approach then directly gives r in terms of the Green's function. This second approach is 
more powerful since it can be used for arbitrary harmonic structures where it is not possible 
to think of simple plane wave scattering states. This approach tells us that we need to 
construct scattering states by evolving the eigenmodes of the two isolated reservoirs. Indeed 
this is what the NEGF approach does in effect and our explicit calculations for a simple 
but representative model clarifies the picture. Our exact calcultions also illustrate some of 
the subtle points involved, such as the correct computation of the self-energies for 
inhomogeneous chains, appropriate normalizations of normal-modes and the choice of initial 
states. 
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